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The data of the HEIDELBERG-MOSCOW double beta decay experiment for the 
measuring period August 1990 - May 2000 (54.9813 kgy or 723.44 molyears), pub- 
lished recently, are analyzed using the potential of the Bayesian method for low 
counting rates. First evidence for neutrinoless double beta decay is observed 
giving first evidence for lepton number violation. The evidence for this decay 
mode is 97% (2.2(t) with the Bayesian method, and 99.8% c.l. (S.lo-) with the 
method recommended by the Particle Data Group. The half-life of the process 
is found with the Bayesian method to be T°J^ = (0.8 - 18.3) X 10^5 y (95% c.l.) 

with a best value of 1.5 X 10^^ y. The deduced-i value of the effective neutrino 
mass is, with the nuclear matrix elements from El, (m) = (0.11 - 0.56) eV (95% 
c.l.), with a best value of 0.39 eV. Uncertainties in the nuclear matrix elements 
may widen the range given for the effective neutrino mass by at most a factor 
2. Our observation which at the same time means evidence that the neutrino 
is a Majorana particle, will be of fundamental importance for neutrino physics. 
PACS. 14.69.Pq Neutrino mass and mixing - 23.40.Bw Weak-interaction and lep- 
ton (including neutrino) aspects - 23.40.-s Beta decay; double beta decay; electron 
and muon capture. 



The neutrino oscillation interpretation of the atmospheric and solar neu- 
trino data, deliver a strong indication for a non-vanishing neutrino mass. 
While such kind of experiments yields information on the difference of squared 
neutrino mass eigenvalues and on mixing angles, the absolute scale of the neu- 
trino mass is still unknown. Information inm double beta decay experiments 
is indispensable to solve these questions Eillld. Another important problem is 
that of the fundarnfiiital character of the neutrino, whether it is a Dirac or a 
Majorana particle Hfl. Neutrinoless double beta decay could answer also this 
question. Perhaps the main question, which can be investigated by double 
beta decay with high sensitivity, is that of lepton number conservation or 
non-conservation. 

Double beta decay, the rarest known nuclear decay process, can occur in 
different modes: 
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21//3/3 - decay : A{Z,N) ^ A{Z + 2,N-2) + 2e- + 2ve (1) 

01//3/3 - decay : a\z,N) ^ A{Z + 2,N -2) + 2e- (2) 

0j/(2)x/9/3 - decay : A{Z,N) ^ a{z + 2, N ^i) + 2e- + {2)x (3) 

While the two-neutrino mode (1) is aUowed by the Standard Model of 
particle physics, the neutrinoless mode {0i>P(3 ) (2) requires violation of lepton 
number (AL=2). This mode is possible only, if the neutrino is^a Majorana 
carticle, i.e. the neutrino is its owp._aiitiparticle (E. Majorana. 1_Q. Racah 
u, for subsequent works we refer to ErQu, for some reviews see Eljj'E2Hi3) . First 
calculations f>f Qv[3(3 decay based on the Majorana theory have been done by 
W.H. Furry □. 

Neutrinoless double beta decay can not only probe a Majorana neutrino 
mass, but various new physics scenarios beyond the Standard Model, such as 
R-parity violating supersymmetric models, R-parity conserving SUSY models, 
leptpmiatks, violation of Lorcntz-invariance, and compositeness (for a review 
see Ej'EJ'O) . Any theory containing lepton number violating interactions can 
in principle lead to this process allowing to obtain information on the specific 
underlying theory. The experimental signature of the neutrinoless mode is a 
peak at the Q- value of the decay. 

The HEIDELBERG-MOSCOW double beta decay experiment in the 
Gran Sasso Underground Laboratory searches for double beta decay of 
"^•^Ge — >'^'^ Se + 2 + {2v) since 1990. It is the most sensitive double 
beta experiment since almost eight years now. The experiment operates five 
enriched (to 86%) high-purity ^^Ge detectors, with a total mass of 11.5 kg, 
the active mass of 10.96 kg being equivalent to a source strength of 125.5 mol 
^^Ge nuclei. This is the largest source strength in use. 

The high energy resolution of the Ge detectors assures that there is no 
background for a Qufifi line from the two-neutrino double beta decay in this 
experiment. 

The unique feature of neutrinoless double beta decay is that a measured 
half-life allows to deduce information on the effective Majorana neutrino mass 
(m), which is a superposition of neutrino mass eigenstates. The half- life-is 
given, when ignoring contributions from right-handed weak currents, by cHj. 

[7^V2(0+-0+)]-i = G™„,^ (1) 

(m) = |m(i) I -f e'*^ |mg) | + e'*^ \mf} \ , (2) 

Here jmie'* | exp (i0i) (i = 1,2,3) are the contributions to (to) from individ- 
ual mass eigenstates, with denoting relative Majorana phases connected 
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with CP violation and, Cmm denotes a nuclear matrix element, which can be 
calculated. 

The effective mass is closely related to the parameters of neutrino oscil- 
lation experiments. The importance of (to) for solving the problem of the 
structure of the neutrino mixing matrix and in particular to fix the absolute 
scale of the neutrino mass spectrum which cannot be fised hy v - oscillation 
experiments alone, has been discussed in detail in e.g. O. 

In this paper, we present a new, refined analysis of the data obtained in 
the HEIDELBERG-MOSCOW experiment during the period August 1990 - 
May 2000 Ej. The analysis concentrates on the neulrkipless decay mode which 
is the one relevant for particle physics (see, e.g. MM). A descriptimi of the 
IDELBERG-MOSCOW experiment has been given recently in Ej, and in 

Fig. |l] shows the combined spectrum of the five enriched detectors ob- 
tained over the period August 1990 - May 2000, with a statistical significance 
of 54.981 kgy (723.44 molyears), around the Q-value of double beta decay. 
The latter has been determined to be Q^^ = 2039.006(50) keV in a recent 
precision experiment EZl. Fig. ^ shows the spectrum obtained with detectors 
Nr. 1,2,3,5 over the period August 1990 - May 2000, with a significance of 
46.502 kgy. Fig. || shows the spectrum of single site events (SSE) obtained 
for detectors 2,3,5 in the period November 1995 - May 2000, under the re- 
striction that the signal simultanously fulfills the criteria, of all three methods 
of pulse shape analysis we have recently developed llsM and with which we 
operate all detectors except detector 1 (significance 28.053 kgy). Double beta 
events are single site events confined to a few mm region in the detector, 
corresponding to the track length of the emitted electrons, in contrast to e.g. 
multiple scattering 7 - events. From simulation we expect that about 5% of 
the double beta single site events should bejSeen also as MSE. This is caused 
by bremsstrahlung of the emitted electrons E3. 

All the spectra are obtained after rejecting coincidence events between dif- 
ferent Ge detectors and events coincident with activation of the muon shield. 
The spectra which are taken in bins of 0.36 keV are shown in Figs. |l|,^,|| in 
1 keV bins, which explains the broken number in the ordinate. We do the 
analysis of the measured spectra with (Fig. |l|) and without the data of de- 
tector 4 (Figs. ||,||) since the latter does not have a muon shield and has the 
weakest energy resolution. We ignore for each detector the first 200 days of 
operation, corresponding to about three half-lives of ^^Co (ri/2 = 77.27 days), 
to allow for some decay of short-lived radioactive impurities. 

The background rate in the energy range 2000 - 2080 keV is found to 
be (0.17 ± 0.01)events/kgykeV {without pulse shape analysis) considering 
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Figure 1. Sum spectrum of the ^^Ge detectors Nr. 1,2,3,4,5 over the period August 1990 
to May 2000, (54.981 kg y) in the energy interval 2000 - 2080 keV, around the Q^^ value 
of double beta decay (Q^^ = 2039.006(50) keV). The curve results from Bayesian inference 
in the way explained in the text. It corresponds to a half-life T^|^2 = (0-80 ~ 35.07) X 10^^ 
(95% C.I.). 
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Figure 2. Sum spectrum of the ^^Ge detectors Nr. 1,2,3,5 over the period August 1990 to 
May 2000, 46.502 kgy. The curve results from Bayesian inference in the way explained in 
the text. It corresponds to a half-life T0|;2=(0.75 - 18.33) X 10^^ y (95% c.l.). 



all data as background. This is the lowest value obtained in such type of 
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Figure 3. Sum spectrum, measured with the detectors Nr. 2,3,5 operated with pulse shape 
analysis in the period November 1995 to May 2000 (28.053 kg y), in the region of interest 
for the Oi^/3/3 - decay. Q»jj^vents identified as single site events (SSE) by all three pulse 
shape analysis methods E3't3 have been accepted. The spectrum has been corrected for the 
efficiency of SSE identification (see text). The curve results from Bayesian inference in the 
way explained in the text. The signal corresponds to a half-life Tjj2 = (0. 88 — 22.38) X 10^^ y 
(90% C.I.). 



experiments. The energy resolution at the Qf3f3 value in the sum spectra is 
extrapolated from the strong lines in the spectrum to be (4.00 ± 0.39) keV for 
the spectra with detector 4, and (3.74 ± 0.42) keV (FWHM) for those without 
detector 4. The energy calibration of the epqppjiment has an uncertainty of 
0.20 keV. For more experimental details see Il3't3. 

We analyse the spectra shown in Figs. |^,|^,^ with the following methods: 

1. Bayesian method, which is used widely at present (see, e.g. Ell). This 
method is particularly suited for low counting rates, where the data follow a 

son distribution, that cannot be approximated by a Gaussian (see, also 

2. Method recommended by the Particle Data Group @. 

The ~ method is not used since it is expected to have some limita- 
tion for application at low counting rates, as also the Maximum Likelihood 
Method. The limitation is due to the Gaussian approximations inherent in 
these methods. 

To be conservative, we first concentrate on the spectra without pulse 
shape analysis. For the evaluation, we consider the raw data of the detectors. 
After that we shall turn to the analysis of the spectrum obtained with pulse 
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shape discrimination. 

Figs. ^jsjJI show the resuh of a Bayesian peak detection procedure. One 
knows that the hncs in the spectrum are Gaussians with a standard deviation 
a = 1.70 (1.59) keV corresponding to 4.0 (3.7) keV FWHM. Given the position 
of a hne, we used Bayes theorem to infer the contents of the hne and the level 
of a constant background. We first describe the procedure summarily and 
then give some mathematical details. 

Bayesian inference yields the joint probability distribution for both pa- 
rameters. The background level was integrated out. This yielded the dis- 
tribution of the line contents. If the distribution has its maximum at zero 
contents, only an upper limit for the contents can be given and the procedure 
does not suggest the existence of a line. If the distribution has its maximum at 
non-zero contents, the existence of a line is suggested and one can define the 
probability K that there is a line with non-zero contents. It is associated with 
an error interval containing with probability K the true value of the contents. 
The Bayesian probability K agrees with the confidence of classical statistics 
in the case of a Gaussian likelihood function. 

For every energy E of the spectra of Figs. [l]j|,^, we have determined the 
probability K that there is a line at E. All the remainder of the spectrum 
was considered to be background in this search. The result is given on the 
left-hand sides of Figs. |,|,| 

We define the Bayesian procedure in some more detail. It starts from the 
distribution p(a;i...a;M|/9, ??) of the count rates x\...xm in the bins 1..M of the 
spectrum - given two parameters p, j]. The distribution p is the product 

A 

vi:x.x...XM\p-,r]) = Yl—i-e'^" (3) 

fe=i 

of Poissonians for the individual bins. The expectation value is the super- 
position 

Afc = p [v fiik) + (l - v) /2(fc)] (4) 

of the form factors fi of the line and f2 of the background; i.e. fi(k) is the 
Gaussian centered at E with the above-mentioned standard deviation and 
f2(k)=f2 is a constant. 
Since 

M 

J2Mk) = 1 for 1^=1,2, (5) 

fc=i 
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Figure 4. Scan for lines in the full spectrum taken from 1990-2000 with detectors Nr. 
1,2,3,4,5, (Fig. 1^), with the Bayesian method. The ordinate is the probability K that a 
line exists at energy E. Left: Energy range 2000 -2080 keV. Right: Energy range of interest 
around Qpp. 



one has 

M 

Hence, p parametrizes the total intensity in the spectrum, and rj is the 
relative intensity in the Gaussian line. 

Prior distributions are given by Jeffreys' rule (§5.35 oftj). The parameter 
p is integrated out and the posterior distribution P(?7|xi...Xj\/) of the relative 
contents rj of the spectral line is obtained. 

The peak detection procedure yields lines at the positions of known wejak 
7-lines from the decay of ^^"^Bi at 2010.7, 2016.7, 2021.8 and 2052.9 keV B. 
In additijCUjv^ line centered at 2039 keV shows up. This is compatible with the 
Q-value llJ^ES of the double beta decay process. We emphasize, that at this 
energy no 7-line is expected according to the compilations in EH. We do not 
find indif;ations for the lines from ^^Co at 2034. 7keV and 2042 keV discussed 
earlier £3 (Figs. 

Bayesian peak detection suggests a line at Q^^ whether or not one includes 
detector Nr. 4 without muon shield (Figs. The line is also suggested in 

Fig. |6| after removal of multiple site events (MSE), see below. 

On the left-hand side of Figs. the background intensity (l-rj) iden- 

tified by the Bayesian procedure is too high because the procedure averages 
the background over all the spectrum (including lines) except for the line it is 
trying to single out. Therefore on the right-hand side of Figs. the peak 

detection procedure is carried out within an energy interval that does not con- 
tain (according to the left-hand side) lines other than the one at Q/3/3. Still the 
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Figure 5. Left: Probability K that a line exists at a given energy in the range of 2000- 
2080 keV derived via Bayesian inference from the spectrum shown in Fig. ^. Right: Result 
of a Bayesian scan for lines as in the left part of this figure, but in the energy range of 
interest around Clpp. 



interval is broad enough (about ± 5 standard deviations of the Gaussian hue) 
for a meaningful analysis. We find the probability K = 96.5% that there is a 
line at 2039.0 keV in the spectrum shown in Fig. ^ This is a confidence level 
of 2.1 a in the usual language. For the spectrum shown in Fig. we find a 
probability for a line at 2039.0 keV of 97.4% (2.2 a). In this case the number 
of events is found to be 1.2 to 29.4 with 95% c.L. It is 7.3 to 22.6 events with 
68.3% c.L. The most probable number of events (best value) is 14.8 events. 

We also applied the method recommended by the Particle Data Group Ej. 
This method (which does not use the information that the line is Gaussian) 




energy [keV] energy [keV] 



Figure 6. Scan for lines in the single site event spectrum taken frcpiiL 1995-2000 with detec- 
tors Nr. 2,3,5, (Fig. H), with the Bayesian method (as in Figs. M,p|). Left: Energy range 
2000 -2080 keV. Right: Energy range of interest around Q/sp. 
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finds a line at 2039 keV on a confidence level of 3.1 a (99.8% c.l.). Such value 
is obtained in a wide range of the assumed width of the signal and the energy 
range of evaluation. 



N 


N Bin 


Energy keV 


N 


N Bin 


Energy keV 


1. 


5653 


2034.66 


6. 


5666 


2039.33 


2. 


5658 


2036.46 


7. 


5669 


2040.41 


3. 


5660 


2037.18 


8. 


5674 


2042.21 


4. 


5664 


2038.61 


9. 


5680 


2044.37 


5. 


5665 


2038.97 









Table 1. Events classified to be single site events (SSE) by all three methods of PSA, in the 
range 2034.1 - 2044.9 keV (it 3 cr range around Q^^) of the spectrum taken with enriched 
detectors Nr. 2,3,5 in the period November 1995 - May 2000 (28.053 kgy). 

From the analysis of the single site events, we find after 28.053 kgy of 
measurement 9 SSE events in the region 2034.1 - 2044.9 keV (± 3cr around 
Q/3/3) (Table |l]). Analysis with the Bayesian method of the single site event 
spectrum (Fig. as described before, shows again evidence for a line at 
the energy of Q^_a (Fig. H). Analyzing the range of 2032 - 2046 keV, we find 
the probability of 96.8% that there is a line at 2039.0 keV. We thus see a 
signal of single site events, as expected for neutrinoless double betpudecay, 
precisely at the value obtained in the precision experiment of tZI. The 
analysis of the line at 2039.0 keV before correction for the efficiency yields 
4.6events (best value) or (0.3 - 8.0) events within 95% c.l. ((2.1 - 6.8) events 
within 68.3% c.l.). Corrected for the efficiency to identify an SSE signal by 
successive application of all three PSA methods, which is 0.55 ± 0.10, we 
obtain a {)v(3(3 signal with 92.4% c.l.. The signal is (3.6 - 12.5) events with 
68.3 % c.l. (best value 8.3 events). 

The PDG method gives a signal at 2039.0keV of 2.8 a (99.4%). The 
analysis given in Fig. |^ (left part) shows partly drastic differences in the 
relative intensities of other identified lines compared to the full spectra, which 
should reflect their different composition of single site and multiple site events. 

We emphasize that wc find in all analyses of our spectra a line at the value 
of Q/3/3- Under the assumption that the signal at Q/3/3 is not produced by a 
background line of presently unknown origin, we can translate the observed 
number of events into half-lives. We give in Table 2 conservatively the values 
obtained with the Bayesian method and not those obtained with the PDG 
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method. Also given in Table jipare the effective neutrino masses {m) deduced 
using the matrix elements of EJ-Ej. 

We derive from the data taken with 46.502 kg y the half-hfe T'/^j = (O-^- 
18.3) X 10^^ y (95% c.L). The analysis of the other data sets, shown in Table 
^, and in particular of the single site events data, which play an important 
role in our conclusion, confirm this result. 



Significan- 
ce [kgy] 


Detectors 


-'-1/2 y 


(m) eV 


Conf. 
level 


54.9813 
54.9813 
54.9813 


1,2,3,4,5 
1,2,3,4,5 
1,2,3,4,5 


(0.80- 35.07) X 10^5 
(1.04- 3.46) X 10^5 
1.61 X 10^5 


(0.08 - 0.54) 
(0.26 - 0.47) 
0.38 


95% c.L 
68% c.l. 
Best Value 


46.502 
46.502 
46.502 


1,2,3,5 
1,2,3,5 
1,2,3,5 


(0.75 - 18.33) X 10^^ 
(0.98 - 3.05) X 10^5 
1.50 X 10^5 


(0.11 - 0.56) 
(0.28 - 0.49) 
0.39 


95% c.l. 
68% c.l. 
Best Value 


28.053 
28.053 
28.053 


2,3,5 SSE 
2,3,5 SSE 
2,3,5 SSE 


(0.88 - 22.38) X 10^^ 
(1.07- 3.69) X 10^5 
1.61 X 10^5 


(0.10 - 0.51) 
(0.25 - 0.47) 
0.38 


90% c.l. 
68% c.l. 
Best Value 



Table 2. Half-life for the neutrinoless decay mode and deduced effective neutrino mass from 
the HEIDELBERG-MOSCOW experiment. 



The result obtained is consistenLovith the limits given earlier by the HEI- 
DELBERG-MOSCOW experiment Ej, considering that the background had 
been determined more conservatively there. It is also consistent with all other 
double beta experiments - which still reach less sensitivity. A second Ge-expe- 
riment, which has stopped operation in 1999 after reaching a significance of 
9 kg y yields (if one believes their method of 'visual inspection' in their data 
analysis) in a conservative analysis a limit of T^^j > 0.55 x 10^^ y (90% 
c.L). The ^^^Te geochemical experiment yields (m^) < 1.1 eV (68% c.L), the 
^^^Te cryogenic experiment yields (m^) < 1.8 eV and iJip-CdW04 experiment 
(to^) < 2.6 eV, all derived with the matrix elements lifEl'L^I to make the results 
comparable to the present value (for references see E3) . 

Concluding we obtain, with more than 95% probability, first evidence for 
the neutrinoless double beta decay mode. 

As a consequence, on the same confidence level, lepton number is not 
conserved. Further the neutrino is a Majorana particle. We conclude from 
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Figure 7. The impact of the evidence obtained for neutrinoless double beta decay in this 
paper (best value of the effective neutrino mass (m) = 0.39 eV, 95% confidence range (0.05 
- 0.84) eV - allowing already for an uncertainty of the nuclear matrix element of a factor 
of ± 50% on possible neutrino mass schemes. The bars denote allowed ranges of (m) 

different neutrino mass scenarios, still allowed by neutrino oscillation experiments (see 
). Hierarchical models are excluded by the new 0!^/3/3 decay result. Also shown are the 
expected sensitivities for the future potential doub|enbeta experiments CUORE, MOON, 
EXO and the 1 ton and 10 ton project of GENIUS M. 



the various analyses given above the effective mass (m) to be (m) = (O.f f - 
0.56) eV (95% c.l.), with best value of 0.39 eV. Allowing conservatively for an 
uncertainty of the nuclear matrix elements of ± 50% ISnv detailed discussions 
of the status of nuclear matrix elements we refer to l23li3^EMJ) this range may 
widen to (m) = (0.05 - 0.84) eV (95% c.l.). 

In this conclusion, it is assumed that contributions to Qv[3l3 dejcaiz_.from 
processes other than the exchange of a Majorana neutrino (see, e.g. cSM) are 
negligible. 

With the limit deduced for the effective neutrino mass, the 
HEIDELBERG-MOSCOW experiment excludes several of the neutrino mass 
scenarios allowed from present neutrino oscillation experiments (see Fig. |^) - 
allowing mainly only for degenerate and partially degenerate mass scenarios 
and an inverse hierarchy - scenario (the latter being, however, strongly 
disfavored by a recent analysis of SN1987A. In particular hierarchical mass 
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schemes are excluded at the present level of accuracy. 

Assuming the degenerate scenapios to be realized in nature we fix - ac- 
cording to the formulae derived in til - the common mass eigenvalue of the 
degenerate neutrinos to m = (0.05 - 3.4) eV. Part of the upper range is already 

g eluded by tritium experiments, which give a limit of m < 2.2 eV (95% c.l.) 
, The full range can only partly (down to ~ 0.5 eV) be checked by future 
tritium decay exnedments, but could be checked by some future (3(3 experi- 
ments (see, e.g. E3'Lil). The deduced best value for the mass is consistent with 
expectations from experimental fj, cy branching limits in models assuming 
the generating mechanism for the neutrino mass to be also responsible for the 
recent indication for as anomalous magnetic moment of the muon l3. It lies in 
a range of interest also for Z-burst models recently discussed aa_explanation 
for super-high energy cosmic ray events beyond the GKZ-cutoffEj. 

The neutrino mass deduced allows neutrinos to still play an important 
role as hot dark matter in the Universe. 

—New approaches and considerably enlarged experiments (as discussed, e.g. 
in Lj) will be required in future to fix the neutrino mass with higher accuracy. 
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